home *** CD-ROM | disk | FTP | other *** search
/ Sprite 1984 - 1993 / Sprite 1984 - 1993.iso / src / lib / c / etc / RCS / ldexp.c,v < prev    next >
Text File  |  1991-12-03  |  3KB  |  117 lines
  1. head     1.2;
  2. branch   ;
  3. access   ;
  4. symbols  sprited:1.2.1;
  5. locks    ; strict;
  6. comment  @ * @;
  7.  
  8.  
  9. 1.2
  10. date     91.12.02.21.51.16;  author kupfer;  state Exp;
  11. branches 1.2.1.1;
  12. next     1.1;
  13.  
  14. 1.1
  15. date     91.12.02.21.50.38;  author kupfer;  state Exp;
  16. branches ;
  17. next     ;
  18.  
  19. 1.2.1.1
  20. date     91.12.02.21.51.35;  author kupfer;  state Exp;
  21. branches ;
  22. next     ;
  23.  
  24.  
  25. desc
  26. @Source for ldexp() library routine.
  27. @
  28.  
  29.  
  30. 1.2
  31. log
  32. @Fix file name in header comment.
  33. @
  34. text
  35. @/*
  36.  * ldexp.c --
  37.  *
  38.  *      Source code for the ldexp library function.
  39.  *
  40.  * Copyright 1988 Regents of the University of California
  41.  * Permission to use, copy, modify, and distribute this
  42.  * software and its documentation for any purpose and without
  43.  * fee is hereby granted, provided that the above copyright
  44.  * notice appear in all copies.  The University of California
  45.  * makes no representations about the suitability of this
  46.  * software for any purpose.  It is provided "as is" without
  47.  * express or implied warranty.
  48.  */
  49.  
  50. #ifndef lint
  51. static char rcsid[] = 
  52.     "$Header: /sprite/src/lib/c/etc/RCS/ldexp.c,v 1.1 91/12/02 21:50:38 kupfer Exp Locker: kupfer $ SPRITE (Berkeley)";
  53. #endif not lint
  54.      
  55. #include <stdio.h>
  56. #include "errno.h"
  57. /* Largest signed long int power of 2 */
  58. #define MAXSHIFT        (8 * sizeof(long) - 2)   
  59.  
  60. #define    MAXFLOAT    1.7e308
  61.  
  62. double ldexp(value, exp)
  63.      double value;
  64.      int    exp; 
  65. {
  66.  
  67.  
  68.     extern double frexp();
  69.     int    old_exp;
  70.  
  71.         if (exp == 0 || value == 0.0) /* nothing to do for zero */
  72.                 return (value);
  73.         (void) frexp(value, &old_exp);
  74.         if (exp > 0) {
  75.                 if (exp + old_exp > 1023) { /* overflow */
  76.                         errno = ERANGE;
  77.                         return (value < 0 ? -MAXFLOAT : MAXFLOAT);
  78.                 }
  79.                 for ( ; exp > MAXSHIFT; exp -= MAXSHIFT)
  80.                         value *= (1L << MAXSHIFT);
  81.                 return (value * (1L << exp));
  82.         }
  83.         if (exp + old_exp < -1023) { /* underflow */
  84.                 errno = ERANGE;
  85.                 return (0.0);
  86.         }
  87.         for ( ; exp < -MAXSHIFT; exp += MAXSHIFT)
  88.                 value *= 1.0/(1L << MAXSHIFT); /* mult faster than div */
  89.         return (value / (1L << -exp));
  90. }
  91. @
  92.  
  93.  
  94. 1.2.1.1
  95. log
  96. @Initial branch for Sprite server.
  97. @
  98. text
  99. @d18 1
  100. a18 1
  101.     "$Header: /sprite/src/lib/c/etc/RCS/ldexp.c,v 1.2 91/12/02 21:51:16 kupfer Exp $ SPRITE (Berkeley)";
  102. @
  103.  
  104.  
  105. 1.1
  106. log
  107. @Initial revision
  108. @
  109. text
  110. @d2 1
  111. a2 1
  112.  * modf.c --
  113. d18 1
  114. a18 1
  115.     "$Header: ldexp.c,v 1.1 88/09/10 16:39:34 mendel Exp $ SPRITE (Berkeley)";
  116. @
  117.